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We investigate the dynamics of a bulk scalar field with various decay channels in the Randall- 
Sundrum infinite braneworld scenario. A bulk scalar field in this scenario has a quasilocalized mode 
which dominates the late-time behavior near the brane. As for this mode, an interesting point is 
the presence of dissipation caused by the escape of the energy in the direction away from the brane, 
even if the bulk scalar field does not have the interaction with the other bulk fields in the bulk and 
fields on the brane. We can interpret that this lost energy is transferred to the dark radiation. We 
show that such an effective 4-dimensional description for a bulk scalar field is valid including the 
various processes of energy dissipation. 



I. INTRODUCTION 



The braneworld scenario has provided a new scheme 
to derive 4-dimensional effective theory from a funda- 
mental higher dimensional one [1]. In particular, the 
second model proposed by Randall and Sundrum (RS2) 
has demonstrated an exciting possibility that the extra 
dimensions may extend infinitely [2]. A lot of work has 
been done on various aspects of this model. Extensive re- 
search based on relativity and cosmology, which is closely 
related to the issue we will discuss in this paper, has also 
been developed [3-6]. 

The background spacetime of the original RS model 
is given by a Z2-symmetric solution of the 5-dimensional 
Einstein equations with a negative cosmological constant. 
However, in the framework of string theory, it is more 
likely that there are various fields other than the gravi- 
ton in the 5-dimensional bulk. Thus, the dynamics of a 
minimally coupled 5-dimensional massive scalar field in 
the braneworld scenario of RS2 type has been studied, 
and it was shown that a bulk scalar field can play the 
role of inflaton field [7-9]. We refer to such a model as 
the bulk inflaton model. 

Prior to the proposal of the bulk inflaton model, 
Dubovsky, Rubakov, and Tinyakov have analyzed the 
dynamics of a bulk scalar field with the mass m on the 
RS2 Minkowski brane background [10]. They have found 
that the bulk scalar field has a quasilocalized mode which 
dominates the late-time behavior near the brane. When 
this mode dominates, the scalar field seen on the brane 
oscillates with the frequency given by to/ \/2 and decays 



gradually. Besides this interesting intrinsic decay, when 
we try to construct a more realistic bulk infiaton model, 
we also need to take account of the dissipation due to the 
interaction with other bulk fields and fields on the brane* . 
We discuss such decay processes in this paper, focusing 
on the effect on the cosmic expansion law realized on the 
brane. Related discussions about the interaction between 
the bulk and the fields on the brane were given in [11]. 
In particular, generation of dark radiation via collision of 
particles on the brane was discussed in [12]. 

This paper is organized as follows. In Sec. II we will 
review the general formulas for the cosmic expansion law 
under the presence of bulk fields. In Sec. Ill we will study 
various decay channels of the bulk scalar field, and show 
how the decay products affect the cosmic expansion law 
on the brane. Section IV is devoted to summary. 



II. COSMIC EXPANSION LAW IN THE 
GENERALIZED RANDALL-SUNDRUM 
BRANEWORLD 

We consider the braneworld scenario of the RS2 type, 
starting with the 5-dimensional Einstein equations 

Rab - ^gabR + J^59ab = K5 i^ab + SabS{r - Tq)] , (1) 

where r is the proper coordinate normal to the brane, and 
the brane is located at r = tq at the fixed point of the Z2 
symmetry. Sab represents the energy-momentum tensor 
composed of the brane tension a and the contribution of 
matter fields confined on the brane Tab- 



*The latter process was investigated in our previous paper [9]. 
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Sab = -O'Qab + Tab ■ 



(2) peff + 4if peff = [-2 (l + ^) Tabu'^n^ - 2He TabU^n^ 



where qab is the metric induced on the brane. We choose 
a := 6A5/K5 = Q/k'^I so that the Randall- Sundrum 
flat brane is realized as the ground state. Here i is the 
bulk curvature radius. Tab is the effective 5-dimcnsional 
energy-momentum tensor of the bulk fields. Here we do 
not specify the form of Tab- 

It is now well known that the 4-dimensional effective 
Einstein equation is derived from Eq.(l) as [3] 



(4) g = k2 c <T-M<UB), (E). 



(3) 



r=ro 

(8) 



The above Eqs.(5), (6), and (8) completely determine 
the cosmic expansion law once we know the evolution of 
the energy-momentum tensor of the 5-dimensional fields 
evaluated on the brane. 



III. DISSIPATION OF THE BULK SCALAR 
FIELD 



where k\ = K%/i and 



'ill/ 
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Tabq^^q". + qy.. ( T^^n^n" - -T\ 



r=ro 



(E) ^ _±_ 



^a 



^a, b d 

bcdn n q^q^ 



(4) 



where is the unit vector pointing outward and normal 
to the brane. Here T^,y is the energy-momentum tensor 
related to the 5-dimensional Weyl tensor ^^^Cabcd- 

To discuss the braneworld cosmology, we assume the 
spatial homogeneity and isotropy for and for each 
component of the energy- momentum tensor. We define 
the energy density of the 4-dimensional ordinary matter 
fields by p := u^u'^t^^, where m'^ is the timelike unit 
vector associated with comoving observers on the brane. 
Similarly, we define p^'' := u^u'^rj:^} for each component. 
Then, Eq. (3) gives the Friedmann equation on the brane: 



F2 = ^ 



p + p(B)^^iE) 



(5) 



where H is the Hubble parameter on the brane. Here we 
have used the fact p^'"^ = p^/2a. 

Next wc consider the equations for the evolution of the 
energy density. From the 5-dimensional conservation law 
integrated over infinitesimally thin layer near the brane, 
the evolution equation for p is obtained as 



p + 4Hp + Ht"^ = 2TabU''n''\r=ro 



(6) 



where an overdot means differentiation with respect to 

the cosmic proper time t. 

On the other hand, 4-dimensional Bianchi identity im- 
plies 



(7) 



where we have defined := p + p'-'^'> + p(^) + p^^^ and 

rj^"^^ in the same way. Then, using Eq. (6) and the facts 

that r(")",, = (p/a)(T^^ + 2p), r^^)''^ - 2i?T„bn'^n^ and 
^(-E)^i^ = 0, we find that the effective energy density of 
the bulk field peS := 



(B) _^ p{E) satisfies 



As a phenomenological model of a bulk scalar field with 
dissipation, we consider the field equation given by 



N{r) 



(9) 



where rn? is the mass squared of a bulk scalar field and 
N{r) is the lapse function. Here represents the rate 
of the decay on the brane, while Vb{r) is a function spec- 
ifying the rate of energy release to the bulk. We call this 
(j) field the bulk inflaton, although our main interest is 
not in the period of infiation but in the period of reheat- 
ing after inflation. For a moment, we do not specify the 
hmctional form of Tb{r). In this setting, we study the 
cosmic expansion law by analyzing the set of equations 
(5), (6), and (8). 



A. Dark radiation from escaping bulk inflaton 

First, we consider the simplest case in which there is 
no bulk field other than the bulk inflaton field (j) that has 
no direct interaction to the fields on the brane. Namely, 



Tab = T'ab and F^ = F^ = 0, where 



T. 



ab 



',aY,b - Qab [ -j^g'"'4>,c(t>,d + -j^^'f 



(10) 



is the energy-momentum tensor of the bulk inflaton field. 
Examining this simplest case is instructive to figure out 
an important difference between a bulk field and a field 
on the brane. 

According to the analyses of the dynamics of a massive 
bulk scalar field on the RS2 flat brane background [10] 
and on the de Sitter brane background [7], there is a 
quasilocalized mode which dominates the late-time be- 
havior. This quasilocalized mode corresponds to a pair 
of the flrst poles appearing in the retarded Green func- 
tion in the frequency domain. To summarize the re- 
sults obtained in the above references, (j) behaves as 

^ g[±ia;-(3/f/2)-(r/2)]t ^j^j^ 



UJ 



m; 



eff 



Q 



"leff, 



2 



r = V 

■1- -L 



under the conditions 



■Im 



1 iu! 
2~ll 



m£ and Hi «; 1. 



(11) 



(12) 



Here m^g. = and i^C^) = r'(^)/r(-2^) is the poly- 

Gamma function. The bulk inflaton oscillates with the 
frequency u) and decays intrinsically even if we do not 
introduce dissipation effects due to the couphng to other 
fields. From the 5-dimensional point of view, this de- 
cay of the amplitude of oscillation occurs as a result of 
the escape of the energy in the direction away from the 
brane. We note that the above expressions are valid only 
when mlg/H'^ > 9/4. Here our discussion is therefore 
restricted to the reheating era after inflation. 

Let us consider the effective energy density of the bulk 
inflaton field from the viewpoint of 4-dimensional ob- 
servers. In the present case, a Z2-symmetric bound- 
ary condition implies ^^a'^" = 0. Therefore we have 

T^f\°'n^ = 0, and then the energy conservation equa- 
tion (6) reduces to the same equation as in the standard 
4-dimensional cosmology: 



p + 4Hp + Ht^^ = 0. 



(13) 



Furthermore, from Eq. (8), we find that Peff = p^^^+p^'^^ 
satisfies 



Peff + 4iJpeff = -Hi 



m 



(14) 



This equation tells us that the evolution of peff depends 
on the evolution of (j) on the brane. Therefore, in order 
to determine the cosmic expansion law, we need to know 
the evolution of the bulk scalar field. 

Next, for comparison we propose a 4-dimensional 
model composed of a scalar field $ and the dark radia- 
tion in addition to the ordinary matter fields. We should 
stress that in the 4-dimensional model the dark radiation 
simply means a component of energy density decoupled 
from the ordinary matter fields. As for the ordinary mat- 
ter fields, we use the same notations as in the case of the 
5-dimensional model. Then the conservation law for the 
ordinary matter fields is identical to Eq. (13). 

We define a scalar field $ which reproduces the late- 
time behavior of the value of (j) on the brane besides an 
overall normalization factor. We assume that $ has mass 
rrieff and decays only into the dark radiation with the 
decay width resc/2. The equation of motion for the ho- 
mogeneous part of the scalar field is given by 



$ + (3/7 + Fesc)^ + m2ff$ = 0. 



(15) 



Using Fcsc ^ JTi, which follows from m£ <C 1, it is easy to 
see that $ also behaves as oc e[±»'^-(3^/2)-(resc/2)]t ^^^^^ 
the same w in Eq. (11). 



The Hubble equation is given by = {k1/3){p+ pcff), 
with pcs := p'-*-' + p^^\ where p*^*^ and p'--^^ are the 
energy density of the $ field and the dark radiation, re- 
spectively. Since the scalar field $ and the dark radiation 
are decoupled from the ordinary matter fields by assump- 
tion, they separately satisfy the conservation law. The 
conservation law for i> and the dark radiation reduces to 

+ 4Hp,s = -Ht^^^^^ = -H ($2 _ 2m^ff$2) . (16) 



Then, with simple identification 

$ = V^(/), (17) 

it is easy to see that the 4-dimensional and the 5- 

dimensional models presented above realize the same cos- 
mic expansion in the law energy regime (12), in which the 
contribution from p^'^^ is neglected. 

From the viewpoint of the 4-dimensional model, we can 
easily write down the evolution of the energy density of 
$ field as 



(18) 



Then, with the aid of Eq. (16), we have 

p(^)+4ifp(^)=re3c^^ (19) 

Thus, we can naturally interpret that the $ field decays 
into the dark radiation. Therefore, from the correspon- 
dence between the 4-dimensional and the 5-dimensional 
models, we are led to the interpretation that the dark 
radiation is generated as a result of this decay process in 
the 5-dimcnsional model. 

Usually, when we consider decay process, it leaves some 
decay products. However, the present decay process of 
the bulk inflaton apparently does not have corresponding 
decay products, since we do not consider the interaction 
with the other fields. Then, one may think that once the 
bulk field decays in this manner, its effective energy den- 
sity is simply lost without being transferred into another 
form. As was expected from the energy conservation law, 
this naive speculation is incorrect. In fact, the lost en- 
ergy was shown to be transferred into the dark radiation 
component. 

Before closing this section, we mention that p^*' in the 
4-dimensional model does not correspond to p^'^^ = p^^^ 
in the 5-dimensional model. Hence, neither does p*-^-* 
to p(^) since p(*) + p(^) = peff = Peff = p('^) + p(^). 
The fact that p^*^ ^ p^"^^ will be easily seen by ex- 
plicitly evaluating p^^^ by using Eqs.(4) and (10) as 



(£/4)[<^2 ^ ^2^2] ^ (l/4)$2 ^ (to2^/2)$2_ This 



expression shows that p^*^ = p^"^^ + £^'^/4. 
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B. Decay of the bulk scalar field to other fields 

In this section, we analyze the field equation (9) with- 
out neglecting the dissipation terms on the background 
of the 5-diniensional anti-de Sitter space with a bound- 
ary de Sitter brane. We make use of the fact that the 
late-time behavior is determined by the singularities in 
the retarded Green function in the frequency domain, as 
was done in Ref. [9]. 

The background metric can be written as 



dr^ + e^"" dajfg)) , (20) 



where R{y) = £sinh^^(|y| -|- yo). The brane is located 
at y = 0, and yo is specified by sinhyo = H£. The 
proper cosmological time on the brane t is related to t 
by r = Ht. Setting ^ = R-^/^e-^''/H{y,T), Eq. (9) 
reduces to 



[-dl + dl - W{y) - U{y)dr\ * = , 



(21) 



with 



W{y) = 



15 + Am'^e 
4sinh2(|2/|+yo) 



3V1 + 3 
— 6iy) - -U{y), 



U{y)^TdS{y)+h{y)R{y). 



(22) 

(23) 



According to Ref. [9], the poles of the retarded Green 
function correspond to the solutions of Eq. (21) that sat- 
isfy the outgoing boundary condition. The solutions can 
be found in the form of 



(24) 



Then, Eq. (21) reduces to an eigenvalue equation for Up 
having as the eigenvalue. When Fb — 0, the equa- 
tion for Up can be solved with the outgoing boundary 
condition as 



My) = jfK-i/2[^oH\y\ + yo)]- 



(25) 



where v := \^mPP + 4. The normalization constant J\f 
is to be chosen so as to satisfy dylup]"^ = 1. Here 
we have not yet imposed the boundary condition on the 
brane. The condition of the Z2-symmetry 



^-+2TZ 



J y=0+ 



(26) 



gives the equation that determines the eigenvalue, i.e., 
the location of the poles in the retarded Green function. 
We focus on the eigenmode with the eigenvalue p that 
has the largest imaginary part in order to discuss the 
late-time behavior of the bulk inflaton field. Under the 



conditions (12), Eq. (26) for this eigenmode reduces to a 
quadratic equation as 



p + \ip + 



Td 9 
— + - 



m 

2IP 



0. 



(27) 



In order to illustrate the effect due to the dissipation 
in the bulk, we consider two specific analytically solvable 
examples given by 

h{y)R{y) = TbH-^ , [case (I)] , (28) 

h{y)R{y) = 2rbHR{yf , [case (II)] , (29) 

where Fh is a constant. In case (I) the first term p^ in 
Eq. (27) is replaced with p'^ + {ip + 3/2)r,,i/-i, while 
in case (II) in the same equation is replaced with 
m? — 2{ip + 3/2)F(,i?. Then, in both cases we have the 
same equation to determine the eigenvalue of the largest 
imaginary part, and we find that a pair of such eigenval- 
ues is given by 



1 



.(rd + Tfe) 



± 



(3iJ + rd + r(,)2 



(30) 



These eigenvalues show that the time evolution of the 
bulk inflaton field evaluated on the brane at late times 
is identical to that of a 4-dimensional scalar field satisfy- 
ing the equation obtained by replacing V^^sc in Eq. (15) 
with Ftot = Eosc + Erf + F;,. Here we have added Fosc 
which comes from higher order corrections neglected in 
Eq. (27). This result shows that all three decay chan- 
nels damp the oscillation of the bulk inflaton fleld as is 
expected. 

So far we have studied just two speciflc examples. In 
more general cases, the eigenfunction cannot be con- 
structed analytically, but we can justify the extensive 
use of the above results by perturbatively evaluating the 
shift of this pair of flrst eigenvalues. We consider that 
the unperturbed eigenfunction Up{y) is given by that for 
U{y) = 0. Then applying the standard perturbation the- 
ory, we find that the shift of the eigenvalue p^ due to the 
U{y)-ievm. is given by 



Ap^ = -[ip + 



3 



j dyU{y)\up{y)\'' 

R{y)h{y) , 



dy 



I 



dy 



(31) 



{y + Hir 



where we have used the fact that |wp(y)P behaves as 
R{y)^ ~ [^/{\y\ + H£)]^ near the brane. If we define 
— {ip + 2>/2)TbH~^ by the right hand side of the above 
equation, we recover the same results obtained above in 
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the analytically solvable cases. From the above equation, 
we also find that the dissipation is more efficient near the 
brane due to the effect of the warp factor. This fact im- 
plies that the energy transferred to the other bulk fields 
is also possibly trapped around the brane as in the case 
of the (j) field. Of course, however, if the bulk fields do 
not have degrees of freedom localized near the brane, the 
energy density near the brane will soon vanish. 

Now let us consider the cosmic expansion law in the 
present case. We divide Tab into two parts as 



l-ab 



-^ab 



ab 



(32) 



Here T^f^' is the energy-momentum tensor of the other 
bulk fields, into which the bulk inflaton field (j) may de- 
cay. Then, the formulas for the energy density evolution 
of the matter fields on the brane (6) and that of the bulk 
fields (8) reduce to 



p + 4:Hp + Ht"^ = rd$2 + 2T^*^u°n'' 

Peff + 4iJ/9eff = 



(33) 



(34) 



j/=0 



where we used $ defined by Eq. (17). The cosmic expan- 
sion law depends on the nature of the fields composing 



r^^' . Below we discuss various cases one by one. 





Let us begin with the simplest situation that the bulk 
excitations caused by the decay of the (j) field are not 
trapped around the brane. In this case T^j'^|y=o will 
rapidly diminish, and hence Eqs. (33) and (34) become 
almost identical to those in the case without dissipa- 
tion in the bulk. The main effect of Ff, arises in Ftot, 
and it can be absorbed by the replacement of Fesc with 
Fesc + Ff,. Therefore we can immediately interpret that 
the <I>-field decays into the dark radiation with the decay 
rate Fggc+Fft and into the matter fields on the brane with 
the decay rate F^. If the matter fields on the brane stay 
radiationlike, the ratio of the dark radiation compared to 
the ordinary radiation is given by (Fesc + ^b)l^d- This 
ratio is constrained by the nucleosynthesis bound. If p de- 
creases slower like dust matter or cosmological constant 
for a certain period, the constraint by the nucleosynthesis 
becomes weaker. 



2) T^;^\y=o hut T^l'u-nX^o = 
Next we consider the case that a part of resulting bulk 
excitation is trapped around the brane. Then, T^f, \y=o 
will not diminish rapidly. Even in this case, if there is no 
direct interaction between the fields on the brane and the 
bulk fields other than (j), we can set T^f,'u""n"\y=o = 0, 
since this term represents the energy transfer between 



(*). 



them. Therefore the effect of nonvanishing T^^ \v=n arises 
only in Eq. (34). Here the model qualitatively bifurcates 
into two cases, depending on whether this bulk excitation 
is completely trapped or not. 

In the former case, the energy density of the trapped 
excitation does not decay into anything. We cannot fur- 
ther analyze the evolution of this energy density without 
specifying the details of the model. In the latter case, the 
energy of the trapped excitation eventually escapes to the 
infinity, and as a result T^f^' n°'n''\y=o vanishes in the end. 
As in the case of the (j) field discussed in the preceding 
section, vanishing of T^^^n^n'' |y=o will be interpreted as 
the transfer of the energy stored in the bulk modes into 
the dark radiation. As before, since the details of the 
model are not provided, we cannot analyze the evolution 
of the energy density of this trapped excitation before it 
is transferred into the dark radiation. 

As a concrete example, let us consider the case that 
the decay product in the bulk is composed of the first 
pair of eigenmodes of another bulk scalar field. Then, 
since the effective energy density of this mode p*^*' be- 
haves like dust matter as in the case of the cj) field, this 
field tends to dominate the energy density as a result of 
the cosmic expansion if the fields on the brane behave 
as radiation. Therefore, models with bulk fields of this 
type will produce large amount of dark matter, which is 
possibly transferred to dark radiation, and hence will be 
strongly constrained by the nucleosynthesis bound. 



Finally, we consider the case that there is bulk excita- 
tion trapped near the brane and it has direct interaction 
with the fields on the brane. In such models, most of the 
discussions given in the previous case apply in the same 
way. However, the history of the universe can be quite 
different, if this interaction is sufficiently strong. In this 
case, the trapped excitation can efficiently decay into the 
ordinary matter fields on the brane. Then, the fraction of 
dark component is reduced, and hence such models can 
avoid the conflict with the nucleosynthesis bound. 



IV. SUMMARY 

We investigated braneworld models of Randall- 
Sundrum type with a massive bulk scalar field, which we 
call the bulk inflaton here, taking into account various de- 
cay channels. We focused on how the cosmic expansion 
law realized on the brane is affected by the influence of 
the oscillating bulk inflaton and its decay products. The 
cosmic expansion law depends not only on the energy- 
momentum tensor of the fields localized on the brane 
but also on that of the bulk fields. Hence, in order to 
determine the cosmic expansion, we investigated the dy- 
namics of the bulk inflaton including dissipation to the 
matter fields on the brane and in the bulk by using the 
Green function method developed in Refs. [7,9,10]. 
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It was argued that the effect of the bulk inflaton on 
the cosmic expansion can be mimicked by an effective 4- 
dimensional inflaton as long as H^t,^ ^ I and |m^|i!^ <C I 
in Ref . [7] . Examining the late-time behavior of solutions 
of the bulk scalar field, we have shown that this claim is 
valid even if the bulk inflaton has various decay channels. 
The decay of the bulk inflaton occurs not only due to the 
dissipation to other fields on the brane or in the bulk but 
also due to the escape of energy to the direction away 
from the brane. From the analysis of the cosmic expan- 
sion law, we have shown that the decay product of this 
latter intrinsic dissipation should be interpreted as the 
dark radiation from the viewpoint of 4-dimensional effec- 
tive theory. Further, we stress that this dark radiation 
is not identical to the contribution from the projected 
5-dimensional Weyl tensor. 

As for the dissipation due to interaction with matter 
fields on the brane, it has been already shown in Ref. [9] 
that the energy lost by this decay process can be inter- 
preted as simply transferred to the matter fields on the 
brane, which is analogous to the usual decay process of 
the 4-dimensional inflaton field. The consequence of the 
decay due to interaction with other bulk fields has vast 
varieties. If the excitation generated by this decay pro- 
cess is not trapped near the brane at all, the effect is 
equivalent to the case of intrinsic dissipation. However, 
we also found that the dissipation to the bulk is more 
efficient near the brane. Hence, it is likely that a part of 
the lost energy is transferred to a mode localized near the 
brane. In this case, the decay product rather behaves as 
a 4-dimensional field as in the case of the bulk inflaton. 
The effective equation of state of such excitation mea- 
sured through the cosmic expansion law depends on the 
details of the model setup, but at least we know that the 
effective equation of state is dustlike if the decay product 
consists of another massive scalar field. If the equation 
of state is dustlike, the fraction of the energy density of 
this excitation continues to increase during the radiation 
dominant phase. In the case that this excitation does not 
decay into the fields on the brane, it behaves as "dark 
matter" . If this excitation is not completely trapped in 
the neighborhood of the brane, it finally decays into the 
dark radiation in the same way as we have seen for the 
bulk inflaton. This further decay transfers the energy 
from this "dark matter" to the dark radiation, but the 
dark component anyway in total dominates the energy 
density unless the decay rate of this excitation is suffi- 
ciently high or its initial amplitude is extremely small. 
Hence, such kind of models will be severely constrained 
by the requirement of the big-bang nucleosynthesis. 

If the excitation in the bulk generated by the bulk in- 
flaton decay has interaction with the fields on the brane, 
the energy c;an be transferred to the ordinary radiation 
fluid. If such decay process is efficient enough compared 
to the other decay channels, the final fraction of the dark 
component will not be significantly large. Such kind of 
models will be viable. 

In conclusion, we found that the bulk inflaton model 



can reproduce the standard cosmology, but the branch- 
ing ratio of the bulk inflaton decay and the nature of the 
decay products arc constrained. The bulk inflaton model 
might leave observable signatures of the braneworld, and 
it may also provide a key to understanding the origin of 
dark matter. 
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